Emergence of One-Dimensional Physics in the Distorted Shastry-Sutherland Lattice 
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Motivated by the on-going investigation of SrCu2(BOa)2 under pressure, we study a variant of the two- 
dimensional Shastry-Sutherland (SS) spin-1/2 model with two types of dimers. Combined with the frustration 
of the SS model, this modification induces, in a large parameter range, a dimensional reduction at low energies, 
with nearly decoupled effective 5 = 1 Haldane chains forming along one of the diagonals of the lattice. We 
also present evidence that the intermediate plaquette solid phase of the undistorted SS model remains stable in 
a finite region of the phase diagram. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.10.Kt, 75.30.Kz 



Introduction. — An exciting route for exploring new phases 
in frustrated magnets consists in applying hydrostatic pres- 
sure so as to considerably modify the dominant magnetic in- 
teractions and thus alter the relevant low-energy degrees of 
freedom. The quasi-2D compound SrCu2(BC>3)2 [1] offers 
such a possibility since it is in the close vicinity of a quan- 
tum phase transition from an exact tensor product of singlets 
toward a plaquette solid[2, 3]. Indeed, Nuclear Magnetic Res- 
onance (NMR)[4, 5], Inelastic Neutron Scattering (INS) [6], 
as well as earlier susceptibility measurements[7], show strong 
evidence that SrCu2(BC>3)2 undergoes at least one phase tran- 
sition under high hydrostatic pressure, the nature of which is 
still not fully understood. What is clear however is that the 
four-fold (C4) symmetry around the void plaquettes is quickly 
lost under pressure, so that nearest-neighbor dimers pointing 
in different directions become inequivalent (see Fig.l). This 
naturally leads to an extension of the Shastry-Sutherland (SS) 
model with two different diagonal bonds Ji and J2 described 
by the spin S — 1/2 Heisenberg Hamiltonian 
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where J' denotes the nearest-neighbor exchange (see Fig. 1). 
Since all nearest-neighbor bonds J' are equal, the product of 
singlets on J\ and J2 bonds is always an eigenstate, as it is for 
the undistorted SS model. Besides, although the C4 symmetry 
is lost, both the sets of J\ and J2 bonds form a square lattice. 
It thus seems natural to expect that the phase diagram will 
consist of three phases, as for the regular SS model: the exact 
singlet phase for small J', the Neel phase for large J', and an 
intermediate plaquette phase with 2D character[2, 3], 

In this paper, using a variety of analytical and numerical 
techniques, we show that making the diagonal bonds inequiv- 
alent actually triggers an abrupt dimensional reduction into a 
new, effectively ID Haldane phase in a wide parameter range 
(see Fig. 2). In this phase the effective low-energy physics 
can be described in terms of nearly decoupled, fully frustrated 
2-leg ladders (see Fig. 1) running along one of the two di- 
agonals of the square lattice. Remarkably enough, the spins 
on the ,/2-bonds form triplets in this phase, leading to a com- 



pletely different picture of the low-lying magnetic excitations 
as compared to the exact singlet phase. 

The emergence of this new singlet phase is best understood 
by starting from the limit J2/J1 = J'/Ji = 0. The model 
then reduces to a set of isolated dimers on the J\ -bonds and 
isolated free spins on the remaining sites (see Fig. 1), result- 
ing into a macroscopically degenerate ground state manifold. 
The degeneracy is lifted by switching on the inter-dimer cou- 
pling J' which, to leading (second) order, a priori couples a 
given spin to 12 neighbors (see Fig. 1, where the 12 neighbors 
of site are numbered from 1 to 12). Each second-order pro- 
cess involves an intermediate state with a J% dimer in a triplet 
state, and the resulting coupling is ferromagnetic and of am- 
plitude — J' 2 /2Ji if the J' bonds are connected to the same 
site of the J\ dimer, while it is antiferromagnetic and of am- 
plitude J' 2 /2Ji if they are connected to opposite sites. When 
processes with different signs couple two spins, they cancel 
out due to quantum interferences. As a result, the couplings 
between site and sites 6 to 12 vanish, and site is only 




FIG. 1: (color online) The distorted SS model (thick and thin gray 
lines for Ji, J 2, and dashed for J') of this work, together with the 
effective quasi- ID model (red solid and dashed thick lines) emerging 
in the limit Ji S> J2, J' which gives rise to the Haldane phase. 
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FIG. 2: (color online) The phase diagram of the present model 
obtained by exact diagonalizations (symbols with lines). The rt 
lines enclose the stability region of the PRVB phase as obtained fro 
the quadrumer boson operator method (see text). The dashed lin 
LI and LI' denote the boundary between the exact dimer and tl 
Haldane phase as obtained from perturbation theory up to second 
and third orders, while the dashed line L2 is the line along which we 
examine the low-E spectrum shown in Fig. 5 below. 



coupled to its neighbors on the ladder shown in Fig. 1. The 
couplings between site and sites 2, 3, 4 and 5 are all anti- 
ferromagnetic and of amplitude J' 2 / 2 J\, while the coupling 
between and 1 is ferromagnetic and equal to — J' 2 / J\ . 

Including both J 2 and J', and up to second order in J' / J\, 
our 2D model thus reduces to an infinite set of 2-leg, frus- 
trated spin ladders directed along the J\ diagonal bonds. The 
effective interactions within the ladders are depicted in Fig. 1 
on top of the original model. They consist of the rung cou- 
pling Jr = J2 — J' 2 /Ji, the leg coupling Jl = J' 2 /2J\, 
and a frustrating coupling J x = Jl- This model has been 
studied quite extensively [8, 9]. It has the special property that 
the total spin on each rung is a good quantum number, which 
allows to determine a large number of eigenstates exactly[8]. 
In particular, in the limit Jr 3> Jl (or J 2 3> J' 2 / Ji) the low- 
est energy is obtained by minimizing the number of triplets 
in the rungs [10]. The resulting ground state of the ladder is 
the product of rung singlets, which together with the strong 
Ji singlets correspond to the exact orthogonal-dimer phase of 
the SS model. On the contrary, if Jr is not too large (and 
a fortiori if it is negative), fluctuations between neighboring 
triplets dominate and the lowest energy is obtained by plac- 
ing a triplet on each rung. The ladder then behaves like a 
spin-1 chain, and hence the ground state corresponds to the 
Haldane gapped phase[l 1]. The actual transition between the 
exact dimer and the Haldane phase is known[8] to take place 
at Jr/Jl — 1.4 and it is of first order. Rewritten in terms 
of the original couplings of our model, this boundary corre- 
sponds to the line J 2 = 1.7 J 72 / Ji, which is shown by the 
dashed line LI in Fig. 2. The dashed line L'l, defined by 
J 2 = 1.7J' 2 /Ji + 1.55J' 3 /J 2 , corresponds to the bound- 
ary obtained from perturbation theory up to third order[17]. 
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FIG. 3: (color online) Left panel: Ground state expectation values 
of the four different nearest-neighbor exchange energies Si ■ Sj for 
the 32-site cluster (ei,2 for the Ji,2 diagonal bonds, and e' 12 for J'~ 
bonds on plaquettes which contain Ji,2 bonds). Right panels: The 
corresponding local energy profiles for representative cases inside 
the exact singlet state (above) and the Haldane state (below). Positive 
(negative) exchange values are shown by dashed (solid) lines, whose 
width is proportional to the magnitude of these values. 



Third order processes in J'/Ji also induce some coupling be- 
tween nearest and next-to-nearest neighboring ladders, result- 
ing into a complex network[ 1 8] , but the quasi- 1 D character of 
the (short-range) spin correlations and the finite spin gap seem 
to be a robust feature of this phase, as suggested by exact di- 
agonalizations. 

Exact diagonalizations. — To determine the behavior of the 
present model away from the above perturbative regime and to 
accurately obtain the full phase diagram of the model, we have 
performed an extensive exact diagonalization study (ED) on 
finite-size clusters with periodic boundary conditions. Apart 
from U(l) spin rotations and translational invariance we have 
also exploited the two reflection symmetries (denoted here by 
72-1 and IZ2) along the two diagonals of the square lattice. 

The transition between the exact singlet phase and the Hal- 
dane phase is most simply identified by measuring the ground 
state spin-spin correlations (Sj -Sj) on neighboring bonds, as 
shown in Fig. 3: they are rigourously equal to —3/4 on dimer 
bonds and to zero on J' bonds in the exact dimer phase, and 
jump upon entering the Haldane phase to ferromagnetic and 
nearly equal to 1/4 on the J2 bonds, and to intermediate anti- 
ferromagnetic values on J\ and J' bonds. The resulting first 
order transition lines shown in Fig. 2 are in good agreement 
with the perturbative predictions and show a weak system size 
dependence. In addition, the correlation profiles shown in the 
right panels of Fig. 3 demonstrate how the 2D-lattice decou- 
ples into almost non-interacting 2-leg ladders. 

Next, we examine the behavior of the spin gap by look- 
ing at the evolution of the low-energy excitation spectrum as 
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we move away from the exact dimer into the Haldane phase. 
Figure 4 shows the low-energy spectrum of the 32-site clus- 
ter along the line J' = 0.4Ji in the S z = and 1 sectors. 
All energies are measured from the ground state. In the ex- 
act dimer phase the spin gap scales linearly with J2 and thus 
decreases as we go towards the Haldane phase. This happens 
because a magnetic S = 1 excitation results by promoting 
a rung singlet[19] into a triplet which costs energy J2 (dis- 
regarding the polarization corrections to the spin gap driven 
by J'[l])- By contrast when entering the Haldane phase at 
J2 — 0.365 Ji the gap does not change with J2 any longer. 
This happens because once we pay the energy to form the ef- 
fective spin-1 entities on the rungs the energy scale J 2 does 
not enter their effective dynamics any longer. 

The ID character of the Haldane phase underlies also a 
number of striking spectral features in the singlet sector. In 
particular, Fig. 4 shows a number of singlets (highlighted 
by dashed lines) which group together into almost degenerate 
levels, and all seem to become soft together with the Haldane 
state (red dashed line). These correspond to singlets where 
not all ladders are in the Haldane phase. Specifically, the 32- 
site cluster (enclosed by solid yellow lines in the right panels 
of Fig. 3) decomposes into four 2-leg ladders with two J2- 
bonds each. To form the Haldane phase in a single ladder we 
pay energy 2 J2 to create two triplets (disregarding again the 
polarization energy driven by J'), but we gain an energy Sh 
from their interaction. This gives 4 possible states with en- 
ergy E\h — 2 J2 — Sh whereby one ladder has entered the 
Haldane phase but the remaining three have still singlets in 
their rungs. Similarly, there must exist: 6 states with two lad- 
ders in the Haldane phase and energy E2H — 2-Ei/j; 4 states 
with three ladders in the Haldane phase and E^h — 3Eih\ 
and finally a single state where all ladders are in the Haldane 
phase with E^h — ^E\h- Both the multiplicities and the de- 
pendence on J2 (the slope in the dashed lines of Fig. 4) are in 
agreement with the ED spectra. In addition, we have derived 
the decomposition of each group of singlets into irreducible 
representations of the space group of the 32-site cluster and 
the results also agree with the ED spectra. Finally, the fact 
that E n n — nE\H explains why the above groups of singlets 
seem to become soft together with the Haldane state. 

We note here that the above energy considerations are only 
true to leading order in J'. There are various perturbative cor- 
rections which renormalize the energy of each of the above 
singlets in a different manner. For the 32-site cluster, it is the 
global Haldane state that crosses the exact singlet phase first. 

Upon increasing J'/Ji, the Haldane phase undergoes a sec- 
ond order transition into the Neel phase. The disappearance of 
the Neel order is continuous and corresponds to the vanishing 
of the local moment m obtained by a usual[3] finite-size ex- 
trapolation tun — m + using numerical data for the cor- 
responding spin structure factor at zero momentum. [20] We 
find that for J2 < 0.6 Ji, the Neel boundary does not depend 
on J2, which is related to the absence of this energy scale from 
the effective dynamics of the Haldane phase. On the other 
hand, the bending of the Neel boundary for J2 > 0.6 Ji indi- 
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FIG. 4: (color online) Low-energy spectrum (measured from the GS) 
for the 32-site cluster along the line J' /Ji = 0.4 in the S z = and 
1 sectors. The blue solid line shows the energy of the exact singlet 
phase. The dashed lines highlight the n-Haldane chain singlets dis- 
cussed in the text. The quantum numbers shown in the legends corre- 
spond to the allowed momenta of the 32-site cluster and the parities 
with respect to the two diagonal reflections TZi t 2 . 



N=32, LE-Spectrum along J =2.89J'-1 




FIG. 5: (color online) Low-energy spectrum (measured from the GS) 
of the 32-site cluster along the line L 2 (J2 = 2.89J' - 1) of Fig. 
2 in the Sz = and 1 sectors. The solid (blue) line stands for the 
energy of the exact singlet phase. 



cates the proximity to another phase, different from the Hal- 
dane one. We now show evidence, first by a detailed spectral 
analysis and then by a perturbative mean-field approach, that 
this phase is adiabatically connected to the plaquette phase 
present for J2 = */i[2, 3]. We begin by examining the evo- 
lution of the low-lying spectrum along the line L2 of Fig. 2. 
As shown in Fig. 5 for the 32-site spectrum, a low-energy 
excitation becomes gradually soft above J2/J1 ^0.7 — 0.8. 
This excitation has zero momentum but it has odd parity with 
respect to IZi and IZ2, in contrast to the ground state. If these 
two states become degenerate in the thermodynamic limit then 
the above quantum numbers are exactly the ones expected for 
the Z2 plaquette phase, since the latter does not break transla- 
tional invariance but it does break IZi and 1Z2- 

We also find a large number of singlets below the lowest 
triplet excitation. The majority of these states have a positive 
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curvature with J \, similarly to the exact singlet state shown 
by the solid (blue) line. Hence these states are most proba- 
bly adiabatically connected to singlets above the dimer phase, 
such as the ones we discussed above that must become nearly 
soft at the boundary with the Haldane phase (in addition there 
could also be some bound states of triplons[12]). In agree- 
ment with Lauchli et al. [3], we also find a singlet which lies 
slightly above the lowest two singlets at J% = Ji . The nature 
of this singlet might well be a singlet excitation (e.g. a domain 
wall) of the plaquette phase. 

A complementary confirmation that the plaquette phase sur- 
vives in a small region for J2 7^ J\ can be also given by 
an analytical mean-field approach based on quadrumer boson 
operators[3, 13-15]. Here, one essentially expands around the 
explicitly symmetry broken phase where every second void 
square is in the singlet ground state of an isolated Heisenberg 
plaquette[2]. The ground state and low-lying excitations of 
the original model are then approximated by truncating the 
Hilbert space of each plaquette to its singlet ground state and 
its lowest triplet, which is achieved by the bosonic represen- 
tation of the four spins (j = 1 , . . . 4) 



(2) 



where a = x,y,z, and s^, create respectively the sin- 
glet ground state and the lowest triplet excitation of a single 
Heisenberg plaquette. A hard-core constraint restricts the bo- 
son occupation number of each plaquette s^s + ^ a t' a t a = 1 
so that the bosonic representation obeys the spin algebra. 
The energy of a single plaquette becomes "Hp = EgS^s + 
Et Yl a tata- Following the standard procedure [3, 14, 15], we 
obtain a triplet mode with dispersion (in units of Ji = 1) 



w 2 (k) = J' 2 + -J' (1-2J') cos k x +{J 2 -2J') cos k y 



(3) 



For J' > i Ji and J 2 < 2 J' the minimum of the spectrum is 
at k = and this becomes soft along the line J2 = — 1 + § J', 
which corresponds to the Neel ordering. The first-order transi- 
tion line between the dimer and the plaquette phase is obtained 
by comparing the exact dimer energy with that of the plaque- 
tte state including the zero-point motion from the quadratic 
terms. These lines are depicted as red lines in Fig. 2. The 
stability region is shifted to larger J' /J, as already observed 
for the standard SS model, but the slopes are in good agree- 
ment with ED. Note that the first order boundary between the 
plaquette and the Haldane phase cannot be found with our 
method since we do not have an expression for the ground 
state energy of the Haldane phase. 

Conclusion. — We have introduced a variant of the SS 
model with two types of diagonal couplings as a minimal 
model for the loss of the C4 symmetry in SrCu2(BC>3)2 un- 
der high pressure. We find a novel quasi-lD Haldane phase 



where the low-energy physics can be described in terms of 
fully frustrated decoupled 2-leg ladders which form 5=1 
entities on their rungs. When J2/J1 decreases, the gap to the 
first excitation first decreases linearly, and then levels off af- 
ter crossing the boundary between the exact dimer phase and 
the Haldane phase. Assuming that the main effect of pressure 
is to reduce the ratio J2/J1, this is precisely what has been 
reported by INS [6], as well as by NMR[16] for the nonmag- 
netic sites. However, the presence of magnetic sites possibly 
arranged on every second ladder, as suggested by NMR[5], 
and of a low-lying excitation below the excitation that levels 
off, as revealed by INS [6], goes beyond the prediction of the 
present model and requires further investigation. 
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